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ABSTRACT 

We model the large kinematic data sets for the four Milky Way dwarf spheroidal 
(dSph) satellites: Carina, Fornax, Sculptor and Sextans, recently published by Walker 
et al. The member stars are selected using a reliable dynamical interloper removal 
scheme tested on simulated data. Our member selection is more restrictive than the 
one based on metallicity indicators as it removes not only contamination due to Milky 
Way stars but also the unbound stars from the tidal tails. We model the cleaned data 
sets by adjusting the solutions of the Jeans equations to the profiles of the projected 
velocity dispersion and kurtosis. The data are well reproduced by models where mass 
follows light and the best-fitting stellar orbits are isotropic to weakly tangential, as 
expected from the tidal stirring scenario. The Fornax dwarf, with more than 2400 
member stars, is a dSph galaxy with the most accurately determined mass to date: 
its la error following from the sampling errors of the velocity moments is below 5 
percent. With mass-to-light ratio of 97 solar units, Sextans seems to be the most dark 
matter dominated of the four dSph galaxies. 
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Carina, Fornax, Sculptor and Sextans - galaxies: fundamental parameters - galaxies: 
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1 INTRODUCTION 

Dwarf spheroidal (dSph) galaxies of the Local Group are be- 
lieved to be among the most dark matter dominated objects 

■ in the universe. Their masses have important implications 
for the problem of the missing satellites. The Carina, For- 

, nax, Sculptor and Sextans dSph galaxies, all known for a few 
decades, are classic, well studied examples of this family of 
stellar systems. Although generally believed to be strongly 
dark matter dominated, their actual mass-to-light ratios are 
still very uncertain. While the tidal origin of their large ve- 
locity dispersions is almost certainly ruled out, their forma- 
tion scenarios are still debatable. 

The accuracy of mass determinations of dSph galaxies 
can only be improved by increasing the size of kinematic 
samples used in the modelling. Recently, new measurements 
of radial velocities for the four dwarfs were published by 
Walker, Mateo & Olszewski (2008a), significantly increas- 
ing the number of stars with measured velocities (by an 
order of magnitude in the case of Fornax). In this Letter 
we use these data to constrain the total mass and velocity 
anisotropy of the dwarfs. Our method relies on the modelling 
of the velocity moments, the dispersion and kurtosis of the 
line-of-sight velocity distribution by solutions of the Jeans 
equations (Lokas 2002; Lokas, Mamon & Prada 2005). An 



underlying assumption of such an analysis is that in spite of 
being tidally affected by the gravitational field of the Milky 
Way the dwarfs can still be approximated as self-gravitating, 
virialized systems. An important first step in the analysis is 
to use a reliable procedure for the removal of interlopers: 
the unbound stars from the Milky Way and the tidal tails 
of the dwarfs. We will show that for all the dwarfs the data 
are well reproduced by models where mass traces light and 
the velocity anisotropy parameter is constant with radius. 

All the assumptions adopted here have been verified by 
A-body simulations performed within the framework of the 
tidal stirring scenario (Mayer et al. 2001). Such simulations 
show that dSph galaxies form effectively from initial stellar 
disks embedded in dark matter haloes by tidal interactions. 
The products of such evolution have mass-to-light profiles 
and anisotropy parameters almost constant with radius (Kli- 
mentowski et al. 2007; Mayer et al. 2007). In spite of losing 
considerable fractions of mass by tidal forces, such objects 
are still gravitationally bound and their velocity moments 
can be reliably reproduced by equilibrium models yielding 
realistic mass and anisotropy estimates (Klimentowski et al. 
2007; Lokas et al. 2008). Since tidal tails in the immediate 
vicinity of the dwarfs are most likely oriented along the line 
of sight and are quite dense (Klimentowski et al. 2008), we 
can however expect significant contamination of kinematic 
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Table 1. Adopted parameters of the dwarfs. 





Carina 


Fornax 


Sculptor 


Sextans 


centre (J2000) 


RA=6 h 41 m 37 s 


RA=2 h 40 m 04 s 


RA=l h 00 m 09 s 


RA=10 h 13 m 03 s 




Dec=-50°58'00" 


Dec=-34°31'00" 


Dec=-33°42'30" 


Dec=-1°36'54" 


distance modulus {in — A^)o 


20.03 ± 0.09 


20.7 ± 0.12 


19.54 ± 0.08 


19.67 ± 0.08 


distance D 


101 ± 5 kpc 


138 ± 8 kpc 


79 ± 4 kpc 


86 ± 4 kpc 


apparent magnitude ray 


11.0 ±0.5 


7.4 ±0.3 


8.6 ±0.5 


10.4 ±0.5 


absolute magnitude My 


-9.03 ±0.59 


-13.3 ±0.42 


-10.94 ±0.58 


-9.27 ±0.58 


luminosity Ly 


(3.50 ± 1.35) x 1O 5 L 


(1.79 ±0.69) x 1O 7 L 


(2.03 ±0.79) x 1O 6 L 


(4.37 ± 1.69) x 1O 5 L 


Sersic radius ijg 


5.5 arcmin 


9.9 arcmin 


6.8 arcmin 


15.5 arcmin 


Sersic parameter m 


1.0 


1.0 


1.0 


1.0 


major axis PA 


65° 


41° 


99° 


56° 
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Figure 1. The selection of member stars. Each point shows the 
heliocentric velocity of a star versus its projected distance from 
the centre of the dwarf. Blue dots correspond to the stars accepted 
by our iterative interloper removal scheme in the final iteration. 
Red dots indicate the stars rejected by the procedure. 



samples by unbound stars. Still, such contamination can be 
effectively removed by an iterative procedure based on a 
maximum velocity available to a given star (den Hartog & 
Katgert 1996; Klimentowski et al. 2007; Wojtak et al. 2007). 



2 THE DATA 

Figure [T]presents the kinematic data sets for the four dwarfs 
from Walker et al. (2008a). In each panel the dots correspond 
to heliocentric velocities of the stars as a function of their 
distance from the galaxy centre. The positions of the centres 
and the distances of the dwarfs were adopted from Mateo 
(1998) except for Fornax, where a newer estimate of the 
centre by Walcher et al. (2003), which agrees well with the 
one of Coleman et al. (2005), was used. All these values 
are listed in Table [T] The velocity range shown in Fig. [T] 
corresponds to a cut-off of about 4a with respect to the mean 
velocity of the galaxy, where a is the velocity dispersion in 



the centre, where the contamination is low. Colour-coded in 
the Figure are the results of the interloper removal scheme 
applied to the data. 

The procedure, first proposed by den Hartog & Kat- 
gert (1996) and applied to galaxy clusters, relies on re- 
jecting stars with velocity exceeding the maximum veloc- 
ity available to a star at a given projected distance from 
the centre. The maximum velocity is estimated by assum- 
ing that the star is on a circular orbit with velocity v c ir = 
\GM(r) /r] 1//2 or falling freely in the galaxy's potential with 
velocity Vi n { = \f2v C it . The mass profile needed for the calcu- 
lation of the velocities is estimated in each iteration using the 
standard mass estimator Mvt derived from the virial theo- 
rem (Heisler, Tremaine & Bahcall 1985). The procedure has 
been tested on mock data sets generated from A r -body sim- 
ulations and shown to effectively remove most of unbound 
stars while rejecting very few bound ones (Klimentowski et 
al. 2007; Wojtak et al. 2007). Note that the efficiency of the 
procedure is not restricted to systems where mass follows 
light and it does not impose such models (Klimentowski et 
al. 2007; Sanchez-Conde et al. 2007). 

The results of the application of the procedure to the 
present data are illustrated in Fig. [1] so that the blue dots 
show the stars accepted in the final iteration while the red 
ones indicate those rejected by the procedure. Within the ve- 
locity and distance ranges shown in the Figure we identified 
700, 2445, 1281 and 412 members respectively for Carina, 
Fornax, Sculptor and Sextans. In the following analysis we 
will only use these cleaned (blue) samples. We have verified 
that the procedure does not depend on the initial velocity 
cut-off, i.e. identical final samples are obtained if the initial 
velocity cut-off is 5a or 3a with respect to the galaxy mean. 

Note that our dynamical interloper removal scheme is 
more restrictive (i.e. it rejects more stars) than the method 
of sample selection based on metallicity indicators discussed 
by Walker et al. (2008b). The latter is intended to identify all 
the present as well as former stellar members of the dwarf 
so in principle it rejects only the contamination from the 
Milky Way stars. On the other hand, our scheme rejects 
also the stars that were stripped from the dwarf due to tidal 
interactions, which now populate the tidal tails and are no 
longer good tracers of the dwarf galaxy potential. 

In order to calculate the velocity moments we binned 
the data so that we had 70, 200, 100 and 40 stars per bin 
respectively for Carina, Fornax, Sculptor and Sextans. From 
the data grouped in this way, the profiles of the line-of-sight 
velocity dispersion a(R) and kurtosis n(R) were calculated 
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Figure 2. Velocity dispersion profiles for the member stars in 
the four dSph galaxies. Errors indicate the sampling errors of the 
dispersion. Red lines show the best-fitting profiles from fitting 
velocity dispersion alone, blue ones from fitting velocity dispersion 
and kurtosis. 
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Figure 3. The profiles of the kurtosis-likc variable k for the mem- 
ber stars in the four dSph galaxies. Errors indicate the sampling 
errors. Blue lines show the best-fitting profiles from the joint fit- 
ting of the velocity dispersion and kurtosis. 



using standard unbiased estimators of these moments (see 
Lokas et al. 2005). The dispersion data points a(R) were 
assigned sampling errors of size a / */ 2(n — 1) where n is the 
number of stars per bin. Instead of kurtosis k we use the 
kurtosis- like variable k = (log k) 1 ^ 10 which has a Gaussian 
sampling distribution contrary to the kurtosis k. The values 
of k were obtained with the correction of the standard esti- 
mator of kurtosis K by the bias due to the low number of 
stars per bin so that k = 3K/C where C = 2.75, 2.84, 2.89 
and 2.93 respectively for n = 40, 70, 100 and 200. The sam- 
pling errors of k are 0.018, 0.0146, 0.0124 and 0.00914 for the 
same values of n. These numerical values were estimated by 
Monte Carlo sampling of data from near-Gaussian distribu- 
tion similar to the one described in the Appendix of Lokas 
& Mamon (2003). The profiles of the velocity moments ob- 
tained for the data cleaned of interlopers (blue dots in Fig.[T|) 
are plotted together with their sampling errors in Figs. [2] (ve- 
locity dispersion) and (kurtosis). We can see that they are 
all slowly decreasing with radius so our method of member 
selection leads to significantly different dispersion profiles 
than in the case where no such dynamical scheme is applied 
and the profiles turn out to be flat (Walker et al. 2007). 



3 RESULTS 

The data shown in Figs. [2] and [3] were modelled by adjust- 
ing the projected solutions of the Jeans equations (see Lokas 
2002; Lokas et al. 2005) assuming that the mass distribution 
follows the light and the anisotropy parameter j3 is constant 
with radius. The shape parameters of the distribution of 
light (i.e. the Sersic radius Rs and the Sersic parameter m) 
and the total apparent magnitude were adopted from Irwin 
& Hatzidimitriou (1995) and are listed in Table[T] The Table 



also gives the total luminosity for each galaxy, adjusted to 
the distances adopted here. The listed errors in total lumi- 
nosity include the errors in apparent magnitude as well as 
distance. 

The lines plotted in Figs. [2] and [3] show the best-fitting 
solutions of the Jeans equations. The red lines (for Fornax 
not visible as overplotted by the blue line) correspond to 
the best models when only the dispersion was fitted. The 
blue lines show the results when both velocity dispersion 
and kurtosis profiles are fitted for each dwarf. In each case 
two parameters, the total mass M and the anisotropy pa- 
rameter /3 where adjusted. The best-fitting values of these 
parameters are listed in Table [5] together with la error bars 
following from the sampling errors of velocity moments. The 
Table also gives the mass-to- light ratio with an error includ- 
ing both the inferred uncertainty of the mass and the error 
in total luminosity. 

The obtained constraints on the parameters are further 
illustrated in Fig. [4] which shows the la, 2a and 3a prob- 
ability contours in the M-/3 parameter plane. Again, the 
red contours correspond to the case when only the disper- 
sion profiles were fitted, while the blue ones to the results 
of fitting both velocity moments. The dots of corresponding 
colour mark the best-fitting parameters in both cases. For 
all dwarfs including the kurtosis in the analysis tightens the 
constraints on anisotropy. Except for Fornax, it also shifts a 
little the best-fitting parameters. 

Our total mass estimates are much lower than the virial 
masses found by Walker et al. (2007) from the analysis of 
velocity dispersion using similar data sets, but with different 
member selection and assuming NFW dark matter profiles. 
As already mentioned, their less restrictive member selection 
method leads to flat velocity dispersion profiles which indeed 
can be reproduced by extended dark matter haloes. Our 
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Table 2. Fitted parameters of the dwarfs. 



galaxy fitted /3 A/[1O 7 M ] M/L v [M /L ] X 2 /N 



Carina 


(7 

a + 


K 


-0 04+ - 24 

01+ 019 
u ' ul -0.24 


2.3 ±0.2 
2.3 ±0.2 


67 ±31 
66 ±31 


7.9/8 
12.1/18 


Fornax 


(7 
cr + 


K 


n oo+O.lS 
u -°°-0.19 
-0 32+ 014 


15.7 ±0.7 
15.7 ±0.7 


8.8 ±3.8 
8.8 ±3.8 


6.7/10 
17.7/22 


Sculptor 


(7 
cr + 


K 


-0 09+ 017 
u.uy_ 20 

-U.l»_ .2 


3.1 ±0.2 
3.1 ±0.2 


15.3 ±6.9 

15.4 ±7.0 


14.6/10 
23.5/22 


Sextans 


a 

a + 


K 


22+ 18 

04+ 18 
u - u ^-0.25 


4.0 ±0.6 
4.2 ±0.6 


91 ±49 
97 ±51 


12.7/8 
23.9/18 




p 0- 




1.8 2 2.2 2.4 2.6 2.8 
M [10 7 M ] 



14 15 16 17 

M[10 7 M s ] 




Figure 4. The lcr, 2<r and 3cr probability contours showing con- 
straints on the two fitted parameters: the total mass M and ve- 
locity anisotropy /3 following from the sampling errors of velocity 
moments. Red contours correspond to constraints from fitting the 
dispersion profile alone, while blue ones to the case when both 
moments are fitted. Dots indicate the best-fitting parameters in 
each case. 



results agree however with the masses they estimate for the 
inner parts of the dwarfs. 

In order to illustrate the differences in our approaches 
to interloper rejection we performed a similar modelling for 
Fornax as before but on a sample selected with the tradi- 
tional constant velocity cut-off of ±3<r = ±33 km s _1 with 
respect to the galaxy mean velocity. As expected, the veloc- 
ity dispersion and kurtosis profiles are now flat and con- 
stant M/L models do not provide very good fits. When 
fitting only the dispersion the best-fitting parameters are 
M = (18.6 ±0.8) x 1O 7 M (M/L v = 10.4 ± 4.4M /L ), 
= — 0.77±g;|? with X 2 /N = 21.7/10, when fitting both 
dispersion and kurtosis we find M = (18.4 ± 0.8) x 1O 7 M 



(M/Lv = 1O.3±4.4M /L ), = -0.mZ w with X /N 



58.1/22. The differences with our estimates in Tableware 
not very large, especially when both moments are fitted. 
This is due to the fact that contamination increases the val- 
ues of both dispersion and kurtosis while their dependence 
on anisotropy is different (see the discussion in Lokas et al. 
2008). 

It may seem that the number of stars we rejected as 
interlopers (see Fig. [TJ is large, especially for Fornax. Note 
however, that significant contamination is expected for many 
of the dwarfs since for realistic, eccentric, cosmologically mo- 
tivated orbits their tidal tails are oriented approximately to- 
wards the Milky Way (so towards the observer) for almost 
all positions along the orbit except close to the pericentre, 
where the dwarfs spend only a small fraction of the orbital 
period (Klimentowski et al. 2008). 



4 DISCUSSION 

The new kinematic data sets for the Carina, Fornax, Sculp- 
tor and Sextans allowed us to determine the mass of the 
dwarfs with unprecedented precision. In the case of Fornax, 
with the largest sample, the mass is determined with better 
than 5 percent accuracy. This is due to the very small sam- 
pling errors of velocity moments. The mass-to-light ratios 
are much worse constrained because of rather large errors 
in the determination of luminosity. In spite of these, a di- 
chotomy is seen among the studied dwarfs: while Fornax and 
Sculptor have a mass-to-light of the order of 10 solar units, 
Carina and Sextans have significantly higher values of 66 
and 97. Although these high values still bear an error of 
about 50 percent, they cannot be reconciled with any stellar 
population values thus signifying that these two dwarfs are 
highly dark matter dominated and similar in this respect to 
the Draco dwarf (Lokas et al. 2005). 

Our mass and mass-to-light determination for Carina 
agrees within errors with previous estimates obtained from 
the total velocity dispersion by Munoz et al. (2006) and Ma- 
teo et al. (1993) if we take into account their different as- 
sumptions about the parameters of the luminosity distribu- 
tion. For Fornax the present estimates are somewhat lower 
than the previous ones (Walker et al. 2006a; Klimentowski 
et al. 2007; Lokas, Klimentowski & Wojtak 2007). Lokas et 
al. (2007) found M/L of the order of 11 and isotropic or- 
bits from a similar analysis as presented here but based on 
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~ 200 stars from Walker et al. (2006a) , a number an order of 
magnitude smaller than presently available. The anisotropy 
determined here is weakly tangential; note however that al- 
though the orbits are still quite close to isotropic, isotropy 
is excluded at 3<r level. 

In the case of Sculptor our M/L estimate is an order of 
magnitude lower than the most recent one by Battaglia et 
al. (2008) who find M/L = 158 at the distance of 1.8 kpc, 
but agrees well with earlier estimates by Queloz, Dubath & 
Pasquini (1995) and Armandroff & Da Costa (1986). The 
reason for this large discrepancy lies in the selection of the 
data as well as in models used. Note that the fraction of stars 
discarded by our interloper removal scheme (Fig. [1} is large 
in the case of Sculptor so the data are highly contaminated 
leading to larger velocity dispersion found by Battaglia et 
al. The contamination probably comes from the tidal tails 
of the dwarf: as discussed by Coleman, Da Costa & Bland- 
Hawthorn (2005) the tails are not visible in photometric 
studies so they must be oriented along the line of sight thus 
contributing significantly to the kinematics (Klimentowski 
et al. 2007). In addition, the data sample of Battaglia et 
al. extends to larger distances from the centre of the dwarf 
than the one studied here, making the contamination even 
more probable. This manifests itself in the flat rather than 
decreasing velocity dispersion profile suggesting the presence 
of an extended dark matter halo with a slowly decreasing 
density profile that can indeed lead to the inference of large 
M/L values at large distances. 

The mass-to-light ratio for Sextans estimated from the 
total velocity dispersion by Hargreaves et al. (1994) is 94 
solar units (when translated to our adopted luminosity) in 
perfect agreement with our value of 97. These values are 
significantly higher than those estimated by Suntzeff et al. 
(1993) who find M/L of the order of 30-54. We thus confirm 
that this dwarf is indeed very strongly dark matter domi- 
nated. Our determined mass also agrees with the rough es- 
timate by Kleyna et al. (2004) who found the mass in the 
range 3 x 10 7 — 1.5 x 10 s Mq assuming isotropic orbits in the 
Jeans equation. In agreement with Walker et al. (2006b) we 
find however no evidence for the presence of a cold kinematic 
core in Sextans claimed by Kleyna et al. (2004). 

We conclude that for all four dwarfs the data are satis- 
factorily reproduced by models where mass traces light and 
the anisotropy is close to zero or weakly tangential. Both 
these findings yield support to the formation scenarios based 
on the so-called tidal stirring model (Mayer et al. 2001). The 
model predicts that an extended dark matter halo initially 
present in a dwarf galaxy is strongly stripped due to tidal 
interactions with the Milky Way (Kazantzidis et al. 2004) 
and the final product we see today should have dark mat- 
ter distributed in a similar way as the stars up to distances 
where the tidal tails start to be visible (Klimentowski et al. 
2007; Mayer et al. 2007). The stellar orbits are expected to 
be weakly tangential, e.g. in the simulation of Klimentowski 
et al. (2007) the final stage has mean (5 = —0.13 (see also 
their fig. 4), of the same order as the values estimated here. 
The tidal scenario also predicts that some remnant rota- 
tion could be expected in the dwarfs at present, however, 
contrary to the previous claims e.g. for Carina (Munoz et 
al. 2006) and Sculptor (Battaglia et al. 2008) no significant 
rotation is found in the data studied here. 
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